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A SPACE CRAFT -BORNE GR ADIOMETER 
MISSION ANALYSIS 


P. Argentiero 
R. Garza-Robles 

Goddard Space Flight Center 


INTRODUCTION 

The spacecraft-borne gradiometer is an ideal instrument for globally mapping geopotential 
fine structure. Mean gravity anomalies exhibit considerable orthogonality in gradiometer 
data, and local blocks of gravity anomalies in local blocks of gradiometer data can there- 
fore be successfully estimated (Reference 1). Thus, the data-reduction problem implicit in any 
attempt to obtain a global and detailed gravity-field mapping can be reduced to manageable 
proportions. Gradiometer data are the only known data type which have this property and 
which can be obtained on a global basis. 

To ensure global coverage and sufficient sensitivity, the gradiometer should be mounted on 
a satellite in a polar, circular orbit with an altitude no greater than 300 km. Such orbits 
are difficult to determine. Consequently, to assess the ultimate cost of a gradiometer mis- 
sion, it is necessary to obtain with some precision the orbit-determination requirements of 
the mission. Also, because the instrument output is sensitive to spacecraft orientation, the 
attitude-determination requirements must be specified. In this report, numerical simulations 
are used to determine the effect of orbit and attitude errors on the recovery of gravity 
anomalies from gradiometer data. The instrument simulated is a rotating gradiometer of the 
type under independent development by the Hughes Research Laboratory (Reference 2) 
and the Bell Aerospace Company (Reference 3). The two instruments are electromechanical 
analogs of each other, and their outputs therefore relate to the gravity field in a mathematical- 
ly identical fashion. The Appendix provides a description of the mathematical procedures 
for extracting gravity-field estimates from gradiometer data. 

ORBIT-DETERMINATION REQUIREMENTS FOR A SPACECRAFT-BORNE 
GRADIOMETER MISSION 

Numerical simulation techniques were used to estimate the effect of orbit-determination 
errors on the extraction of gravity-field information from gradiometer data. A general 
description of the simulation procedures follows. 

Attention was focused on an area of the globe between 0 and 60 degrees latitude and 
between 0 and 60 degrees longitude. If a reference field of degree and order 1 2 is used, 
mean gravity anomalies have a standard deviation of 20 milligals, and adjacent anomalies 
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are virtually uncorrelated (Reference 4). Consequently, geopotential fine structure in this 
region was described by four hundred 3-degree gravity anomalies whose values in milligalileos 
(mgal)* were independently chosen from a normal population with a mean of zero and stan- 
dard deviation of 20. A gradiometer was assumed to be mounted on a satellite in a circular, 
polar, 300-km orbit. The sensing plane of the instrument coincided with the orbital plane 
of the satellite. Observations were postulated every 15 seconds, implying a 1-degree latitude 
by 1.2-degree longitude grid of data. The output of the gradiometer is the sum of a com- 
ponent attributable to the reference field and a component attributable to the anomalous 
field as represented by the values of the gravity anomalies. The anomalous portion of the 
gradiometer output was computed by using Stokes’ formula to obtain the contribution to 
the output from each gravity anomaly and summing the results. Instrument noise was 
assumed to be 0.1 etvos units (1 etvos unit = 10' 9 gal/cm), and a random-number generator 
was used to simulate the effects of the noise on the data. 

Computational considerations make it necessary to estimate local blocks of gravity anomalies 
in local blocks of gradiometer data. Reference 3 demonstrates that an efficient estimation 
strategy is achieved if the data block is the same size as, or somewhat smaller than, the 
block of estimated anomalies and if the estimates of anomalies in the outer three layers of 
the block are discarded because of aliasing. This strategy is implemented here. Attention 
is focused on obtaining good estimates of nine centrally located anomalies subtending a 
block 27- to 36-degrees longitude and 24- to 33-degrees latitude. These anomalies must be 
separated from unadjusted anomalies by three layers. Hence, the set of estimated anomalies 
must subtend a block 18- to 45-degrees longitude and 1 5- to 42-degrees latitude. Anomalies 
outside this block are fixed at zero and left unadjusted in the estimation process. Data 
inside a block 21- to 42-degrees longitude and 18- to 39-degrees latitude were processed by 
means of a standard least-squares estimator. Each gravity anomaly outside of the block of 
estimated anomalies was, in effect, fixed at a zero value and left unadjusted. Figure 1 dis- 
plays the results of a simulation when the orbit determination is assumed to be perfect. The 
numbers shown represent the difference in mgal between estimated and true values ot the 
gravity anomalies. The poor estimates of anomalies in the outer layers of the estimated block 
are, of course, attributable to the proximity of unadjusted anomalies. The shaded area repre- 
sents the data block. The nine centrally located anomalies were recovered with an average 
accuracy of 0.025 mm/s 2 (2.5 mgal). 

Orbit-determination errors were simulated by associating with each data point a satellite 
longitude, latitude, and height which were different from the true values. The data which 
are actually processed are assumed to be the difference between the output of the instrument 
and the output predicted by the reference field. Therefore, an orbit-determination error 
leads to the subtraction of the wrong reference field. In the along-track and cross-track 
directions, this effect is negligible; in the radial direction, it is not negligible. At the nomi- 
nal altitude of 300 km, a radial error will cause the reference value of the gradiometer 


^Throughout the text of this document, the measurement unit of milligalileo (mgal) has been converted to a Standard 
International Unit of millimeters/sec 2 (mm/s 2 ). For convenience, the illustrations have not been converted. The simple 
conversion is 1 mgal = 0.01 mm/s 2 . 
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SATELLITE ALTITUDE = 300 km, INSTRUMENT ACCURACY = 0.1 eV 


Figure 1 . Errors in the recovery of 3-degree gravity anomalies 
(in mgal) from gradiometer data. 

output attributable to a 14* degree and order reference field to be biased by 0.001 8 etvos 
units for each meter of error in the radial direction (Reference 4). Therefore, a 50-meter 
radial error causes the reference field to bias the processed measurement by as much as the 
standard deviation of the instrument. The simulations took this effect into account. The 
results were that along-track and cross-track orbit-determination errors lead to very little 
degradation in the ability of a least-squares estimator to extract gravity information from 
gradiometer output. In fact, simulated along-track and cross-track errors of 300 meters had 
no noticeable effect on the recovery of gravity anomalies. Figure 2 shows the sensitivity 
of gravity/anomaly recovery caused by altitude error. The recoverability of gravity anomalies 
was computed by averaging the absolute values of the difference between estimated and 
actual values of the nine centrally located anomalies. Figure 2 shows that altitude errors 
of up to 50 meters cause no serious degradation in gravity anomaly recovery. 

ATT ITU DE-DETERM I NATION REQUIREMENTS FOR A SPACECRAFT-BORNE 
GRADIOMETER MISSION 

A rotating gradiometer senses gradients in a given plane, that is assumed to be perpendicular 
to the spin vector of the satellite. The instrument’s sensitivity to a gravity anomaly will 
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Figure 2. Three-degree gravity anomaly recovery versus altitude 

errror. 

depend on the orientation of the sensing plane. Figure 3 shows the perturbation pattern of 
an 0.01-mm/s 2 (1-mgal) perturbation of a 3-degree anomaly in gradiometer data. Clearly, the 
strongest and most localized observability pattern is obtained when the spin vector has a zero 
elevation. At zero elevation, the perturbation pattern is virtually independent of azimuth 
angle. 

The anomalous part of the gradiometer output is not sensitive to attitude-determination errors. 
Azimuth and elevation errors of as much as 1 0 degrees change the anomalous gradiometer 
output relative to an 0.01-mm/s 2 (1-mgal) perturbation of a 3-degree anomaly by less than 
0.0006 etvos units. A more serious difficulty is that subtraction of the reference field leads to 
a data bias in the case of attitude-determination error. In Reference 4, the gravity-gradient 
tensor of a (14, 14) reference field at altitude 300 km, latitude 37 degrees, longitude 260.5 
degrees is computed. Using this tensor as a reference, the data bias introduced can be com- 
puted by azimuth and elevation error at 300 km when a (14, 14) reference field is used. 

(See Appendix.) A 0.2-degree error in determination of elevation angle leads to a data bias 
of 0.05 etvos units. Results concerning the impact of orbit-determination errors indicate 
that such a bias does not lead to a serious degradation of results and that a 0.2-degree error 
level can be considered a reasonable requirement. The impact of azimuth errors is less severe. 

A 5-degree azimuth error yields a data bias of 0.03 etvos units. Hence, an error level of 5 
degrees is the requirement for azimuth-angle resolution. 
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ELEVATION = 0 



Figure 3. Gradiometer output perturbation caused by 0.01 mm/s 2 (1-mgal) perturba- 
tion of 3-degree gravity anomaly (satellite altitude = 300 km) 


CONCLUSIONS 

Numerical simulations were performed to obtain the orbit- and attitude-determination require- 
ments of a spacecraft-borne gradiometer mission. The simulations were made with the assump- 
tion that the spacecraft was in a circular, 300-kilometer, polar orbit. Observations of a 0.1- 
etvos unit accuracy were assumed every 1 5 seconds. Aliasing effects of distant, unadjusted 
gravity anomalies were included. Results demonstrate that position determination of 
300 meters in the along-track and cross-track directions and 50 meters in the radial direction 
are required for a successful spacecraft-borne gradiometer mission. When these requirements 
are met, the degree gravity anomalies can be recovered with an accuracy of 0.025 mm/s 2 
(2.5 mgal). 

The optimal orientation of the gradiometer sensing plane is achieved when the spin-vector 
elevation is 0 degrees. The attitude-determination requirements are 5-degree spin-vector 
azimuth and 0.2-degree spin vector-elevation. 


Goddard Space Flight Center 

National Aeronautics and Space Administration 
Greenbelt, Maryland June 7, 1976 
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APPENDIX 


MATHEMATICAL PROCEDURES FOR OBTAINING DETAILED GRAVITY 
FIELDS FROM GRADIOMETER DATA 


The following is an outline of the mathematical procedures for using conventional parameter- 
estimation techniques to recover gravity anomalies from the output of a spacecraft-borne 
rotating gradiometer. The next section details the precise relationship between the gradiom- 
eter output and the ambient gravity field. In succeeding sections, a convenient parameteriza- 
tion of the gravity field is displayed, and the sensitivity matrix of this parameterization 
relative to gradiometer data is derived. Finally, an algorithm for combining gradiometer data 
with surface-gravity data and satellite-perturbation data to yield an optimal estimate of the 
gravity field is displayed. 

DERIVATION OF THE OUTPUT SIGNAL OF A ROTATING GRADIOMETER 

Figure A-l represents a rotating gradiometer of the type described in Reference 1 . Assume 
that accelerometers are at positions , P 2 , P 3 , and P 4 and that the system is rotating in the 
plane of the figure with angular velocity co. The output of accelerometer A. at time T is 
<x(T). The four signals are electronically combined to yield an output of the system at time 
Tof 


output = [aj(T) + a 3 (T)] -[a 2 (T) + a 4 (T)l (1) 

Represent the force attributable to gravity in the sensing plane of the instrument as 

F(x.y)-F x l x + F y l y (2) 

The output of accelerometer A. can be expressed in vector cross-produce notation as 

«j(T) = || F(x(T), y(T» ® Pj(x(T), y(T)) || (3) 

For simplicity, assume that P , P 2 , P 3 , and P 4 are unit vectors. Then 

(T) = || F(cos coT, sin coT) ® (cos toT I x + sin coT I y ) || 
a ? (T) = || F(-sin coT, cos coT) ® (-sin coT I x + cos coT I ) || 

(4) 

ck 3 (T) = || F(-cos coT, -sin coT) ® (-cos coT I x + -sin coT I y ) || 
a 4 (T) = || F(sin coT, -cos coT) ® (sin coT I x + -cos coT I y ) || 
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SIGNAL = (A 1 + A 3 ) - (A 2 + A 4 ) 

= 2 (V xx -V yy ) sin2coT-4V xy cos2o;T 
AMP = 2 [(V xx -V yy ) 2 +4V^ y ] ,/2 

Figure A-1 . Output signal of rotating gradiometer. 


From equation 4, 

a 1 (T) = F y (cos coT, sin coT) sin coT - F (cos coT, sin coT) cos coT 
ix y 

a 9 (T) = F (-sin coT, cos co T) cos coT + F (-sin coT, cos coT) sin coT 

2 x y f j ) 

a~(T) - -F Y (-cos coT, -sin coT) sin coT + F (-cos coT, -sin coT) cos coT 

3 x y 

a A (T) = -F (sin coT, -cos coT) cos coT - F (sin coT, -cos coT) sin coT 

4 x y 

Expand F (x, y) in a first-order Taylor series about the origin as 

F x (x, y) = F x (0, 0) + V* x + V* y (6) 

F y (x, y) = F y (0,0)+vyx+yy y 
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where 


dF Y 

v x = — v x = 
x ox y 


OF x OF OF 


Oy 


Oy 


In a conservative force field, V-£-=V-|- Hence, from equations 5 and 6, 


(T) = [F x (0,0) + V x coscoT + V X sin coT] sin coT 

- [F y (0, 0) + V x cos coT + Vy sin coT] cos coT 
« 2 (T) = [F x (0, 0) - V x sin coT + Vy cos coT] cos coT 

+ [Fy(0, 0) -V x sin coT + Vy cos coT] sin coT 

a,(T) = -[F (0, 0) - V x cos coT - V x sin coT] sin coT 
j vv x y 

+ [Fy(0, 0) - Vy cos cjT - V^ sin coT] cos coT 
a 4 (T) = -[F x (0, 0) + V x sin cjT-V x coscoT 1 cos coT 

- [Fy(0, 0) + Vy sin coT - vj cos coT] sin coT 


(7) 


Combining equation 7 with equation 1 yields 


output = 2(V X - V^) sin a coT - 4V X cos a coT (8) 

x y y v J 

The maximum output, which is taken to be the essential measurement of the instrument, 
is 


AMP = 2[(V X - Vp 2 + 4Vy ] 1/2 (9) 

Finally, the differential form on the right side of equation 9 is invariant under rotations in 
the sensing plane of the instrument. This implies, as is necessary, that the gradiometer out- 
put is independent of a particular choice of coordinate set. 

A PARAMETERIZATION OF THE GEOPOTENTIAL FIELD 

To apply standard parameter-estimation techniques to recover the gravity field from gradiom- 
eter data, it is necessary to parameterize the gravity field. However, in this context, the 
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parameterization must be carefully chosen. The global recovery of gravity with any degree 
of detail requires the estimation of many parameters. Unless the parameter set can be de- 
composed into smaller subsets which are separately estimated, the data-reduction problem is 
insurmountable. What is required is a representation of the gravity field, the representation 
of which is perturbed only in a given localized area if a given parameter of the representation 
is perturbed. Because the output of the gradiometer is an in-si tu observation, with this type 
of representation it is possible to estimate local blocks of parameters in local blocks of 
gradiometer data. There are several such parameterizations. The one presented here is con- 
sidered to be the most convenient in terms of optimally combining gradiometer data with 
surface-gravity data and satellite-perturbation data to yield the best estimate of the gravity 
field. 

Represent the geopotential field as 

W = U + T (10) 

where U is a reference geopotential generally defined by a low degree and order spherical 
harmonic expansion and T is the so-called anomalous potential. The anomalous potential 
T at any point above the surface of the Earth can be expressed by means of the discrete 
form of Stokes’ formula as 

T( r , 0, X) = ~ S( r. 0, X, 0|, X[) cos (0j) A0[ AX[ 6gj (0{, X|) 00 

i 

where 5 g. (0.', X.') is a mean gravity anomaly over a block, centered at latitude 0.' and longi- 
tude X.', and referenced to an equipotential surface defined by the nominal field U. The 
expression cos (0.') A0/ AX/ represents the area of the block on which the i* gravity anomaly 
is defined, and S (r, 0, X, 0.', X.') is Stokes’ function given by 

P / p +tcos(0) + D-|\1 

S( r, 0. X, 0 , X ) = t|— + 1 - 3D - t cos (0) ^5 + 3Cn I — J JJ (12) 

where 

0 = cos -1 [sin (0) sin (0') + cos (0) cos (0') cos (X' - X)1 
R 

t = _ 
r 

D = (1 - 2t cos (0) + t 2 )' /2 


and R is the mean radius of the Earth. 
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It will be useful to obtain the gravity-gradient tensor for the anomalous^ field in a local 
topocentric coordinate system. Define an orthogonal coordinate set |l r , 1 0 , where 
1, is directed outward from the Earth, 1 ^ is directed northward, and 1 x is directed eastward. 
Because of the solenoidal and irrotational nature of a conservative force field, the tensor matrix 
is symmetric, and the diagonal elements must add to zero. The tensor assumes the form 


a 2 T 

1 3 2 T 

1 

a 2 T 

3r 2 

r 9r 30 

r 

arax 


1 3 2 T 

1 

a 2 T 


r 2 30 2 


d(j)d\ 


3 2 T 1 3 2 T 

_ 3r 2 r 2 3X 2 _ 

From equation 1 1 , we obtain 


(13) 


a 2 T 

R 

to 

1 

4ir 

a 2 T 

R 

30 2 

47T 

a 2 T 

R 

ax 2 

4n 

a 2 T 

R 

d 09X 

4 tt 

a 2 T 

R 

3r 30 

4tf 

3 2 T 

R 

arax 

4tt 


a 2 s 

aT 2 

a 2 s 


E TT ( r > cos W'O A< ^i AX i 5g i 


E 


a 2 s 


a 2 s 


( r , 0, X, XJ) cos (0j) A0j AXj {$[, X|) 


( r , 0, X, 0[, Xj) cos (0-) A0j AX[ §gj (01, Xj) 


( r, 0, X, 0- , X-) cos (0j) A0j AXj 5g ; (0j, X[) 


(14) 
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An application of the usual chain rule for derivatives yields the following recursion relations 
for the necessary partial derivatives of Stokes’ function 


a = sin 0 sin 0' + cos 0 cos 0' cos (A' - A) 


0 = cos a 


da 

— = cos 0 sin 0 ' - sin 0 cos 0' cos (A' - A) 
d0 

d 2 a 

= _ a 

d0 2 

d 2 a 

= sin 0 cos 0' sin (A' - A) 

dAd0 


da 

cos 0 cos 0 * sin (A' - A) 

— = 

dA 

d 2 a 

- cos 0 cos 0 r cos (A' - A) 

= 

dA 2 

d0 

da 

(1 _ a 2 } -Vi _ 

d0 

d</> 

d 2 0 

/da\ 2 

( — ) ( 1 - a 2 ) 3/2 a + ( 

d0 2 

\ d</>/ 


d 2 0 da da d 2 a 

= - a ( 1 - a 2 )~ 3 ^ 2 — — + ( 1 - a 2 )' 1 / 2 

dAd0 d0 dA dAd0 


d 0 ? 

— - (1 -a 2 )" 1/2 — 
dA dA 


d 2 0 /da\ 2 9 1/9 d2a: 

= — (1 -a 2 r 3/2 a + (l -a 2 r 1/2 

dA 2 VdA / dA 2 


( 15 ) 
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dS 

dip 

a 2 s 

dip 2 


a 2 s 

9 r dtp 

d 2 S 

dr 2 

d 2 S 

dtp 2 

d 2 S 

9X 2 

a 2 s 

dtpd\ 

a 2 s 

dr dtp 

a 2 s 

a rax 


- t 2 sin 

0 

D 3 

r 

2 

- t 2 COS Ip 

— 

.D 3 

+ 1 3 sin 3 ip 


Ld s 


6 

+ - -8-3 
D 


6 

+ - - 8-3 
D 



1 - t cos \p - D 

- 3 2n 

D sin 2 0 

1 - t cos ip - D 

- 3 Cn 

D sin 2 ip 


1 - t cos ip + D 
2 

1 - t cos \p + D 


] 
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The parameterization of the gravity field in terms of gravity anomalies has some advantages. 
The anomalies are, in fact, measurables. Hence, if this parameterization is used to reduce 
gradiometer data, it becomes easy to optimally combine the output of the gradiometer with 
surface-gravity data to yield a best-estimate gravity field. Significant data-reduction advantages 
also exist. Figure A- 2 displays the perturbation pattern of a gravity anomaly in gradiometer 
data. The localized nature of this pattern suggests that if two gravity anomalies are sufficiently 
separated, their perturbation patterns would not overlap, and the anomalies could be esti- 
mated separately without serious aliasing. Local blocks of gravity anomalies can therefore 
be estimated in local blocks of gradiometer data. 

SENSITIVITY MATRIX OF GRAVITY ANOMALIES IN GRADIOMETER DATA 

To apply standard parameter-estimation techniques for recovering gravity anomalies from 
gradiometer data, it is necessary to obtain the so-called sensitivity matrix. The elements 
of this matrix are the partial derivatives of each gradiometer observation with respect to 
each estimated gravity anomaly. 

To derive the sensitivity matrix, choose a coordinate set j I x , I , I z j so that I z is coincident 
to the satellite spin vector. From equations 9 and 10, the amplitude of the output signal 
of the gradiometer can be expressed as 

AMP= 2[(\V XX + T xx - W yy - Tyy) 2 + 4(W xy + T x y) ? I (16) 

where W and T are the second partial derivatives of the reference and anomolous Fields 

XX XX 1 

in the x direction, W and T are the second partial derivatives of the reference and anoma- 
yy yy 

lous fields, and W xy and T xy are the cross-derivatives of the reference and anomalous fields 
in the x and y directions. The derivative of the signal amplitude of a particular point (r, 0, 

A) in space relative to a gravity anomaly centered at (0', A') is 


where 


9AMP 
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3T„ ^ DAMP 3T yy 

a AMP 3T . 
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’ W xx 

3T yy 9g 

3T xy dg 
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v xx yy 



^ T xx 
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u xx yy' xy 
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lvv xx yy' 



aT yy " 
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ii Wyy / vy xy 
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rrw - W ) 2 + 4 W 2 
n xx yy 2 xy 



(17) 


( 18 ) 


16 




Figure A-2. Gradiometer output perturbation attributable to 0.01 
mm/s 2 (1 mgal). 

Assuming that the sensing plane of tile instrument is near its optimal orientation, the above 
derivatives can be computed from the gradient tensor of a (14, 14) reference field provided 
in Reference 2. To five-place accuracy, these derivatives are given as 


3AMP 


3T 


3 AMP 3AMP 

= - 1 . = 0.01 


xx 


3T. 


yy 


9T 


xy 


Hence, to sufficient accuracy, the necessary derivative can be written as 

3AMP _ 3T XX 3T yy 

3g 9g 9g 


(19) 


( 20 ) 
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The derivatives T and T are elements of the gradient tensor of the anomalous field in the 
spacecraft fixed-coordinate set | I x > I y > \ f • Equations 13, 14, and 15 provide the gradient 
tensor of the anomalous field in a local to po centric-coordinate set. Define 

q = gravity-gradient tensor of the anomalous field in a local topocentric- 
coordinate system 

C = gravity-gradient tensor of the anomalous field in the satellite fixed- 
coordinate set 

A = matrix of the linear transformation which maps vectors from the local 
topocentric system to the spacecraft fixed system 

The rules for the tensor calculus assert that 


C = ACA 1 (21) 

If a and / 3 represent the azimuth and elevation respectively of the satellite spin vector, 
then A is 


A = 


cos j 3 -sin /3 sin a 

0 cos a 


cos a 
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-sin (S -cos |3 sin a cos j3 cos a 


( 22 ) 


Equations 13, 14, 21, and 22 yield 

I - / 


3T Rcos0'A0'A\' 


xx 


9g 


47T 


a 2 s 


a 2 s 


cos |3 — - - — sin j3 sin a 

\ 3r 2 r dr d(j> 


I 


i a 2 s \ / i a 2 s 

+ — cos a cos j3 - — cos )3 - 

r 9r 9 X / \r 9r90 r 2 

i a 2 s \ /i a 2 s 

+ cos G! I sin P sin a + — cos jS • 


sin P sin a 


a 2 s 


909A i 


dr dX r' 


90 2 
sin p sin a 


(23) 


9 2 S 


909X 


cos a 


9 2 S 
9 r 2 


I 9 2 S 

COS OL 

r 2 9X 2 


cos a 


9Tyy Rcos0'A0 f AX' 


9g 


47T 


I 9 2 S 

cos a 

r 2 9 0 2 


9 2 S \ 


+ sin a 

r 2 909X 


cos a 


18 



(23 continued) 
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Equations 20 and 23, together with recursion relations 15, provide a computational algorithm 
for obtaining the needed sensitivity matrix. 

MINIMUM VARIANCE ESTIMATOR FOR OBTAINING GRAVITY ANOMALIES 
FROM GRADIOMETER DATA 

'V 

Let 5g be the difference between the true values of the gradiometer observations and 

nominal values as computed from a reference field. Assume that the component of 

6y is the anomalous gradiometer reading obtained when the position of the satellite is 

given in spherical coordinates as (r. , 0., A.). Let g be a fector of numerical values of a 

global set of gravity anomalies. The V th component of g is the numerical value of the 

gravity anomaly centered at latitude 0. and longitude A.'. The functional relationship be- 
^ ^ j 1 

tween 5y and g can be approximated as 

5y = Ag (24) 

where A is the sensitivity matrix, the number of whose rows is the number of observations 
and the number of whose columns is the number of gravity anomalies. The element in the 
\ m row and column is the partial derivative of the i rf7 gradiometer measurement with respect 
to the j tfl gravity anomaly. 

The actual output of the gradiometer, minus the values computed from the reference field, 
provides direct observations 5y of 5y with statistics 

5y = 3y + v s H (v) - 0, H(i^ r ) = Q (25) 

Estimates of mean free-air gravity anomalies obtained from satellite tracking and gravimetry 
measurements are available. Unless this information is correctly factored into the gravity- 
anomaly estimates obtained from' gradiometer data, the resultant estimates are not optimal. 
Consequently, assume the existence of an a priori estimate g of 'g with statistics 

g' = g + a , E(a) = 0, H(aa T ) = Pj (26) 

Under these assumptions, the minimum variance estimator of g is known to be 

'g = (A t Q" 1 A + P" 1 r 1 (A t O" 1 5y + P“ l g) (27) 

And the covariance matrix of this estimator is given by 

E[(t-g)(g-g) T ] = (A T Q- , A + 1>7' r' ( 28 ) 


19 



REFERENCES 


1. Hughes Research Laboratory, “Lunar Orbiting Gravity Gradiometer Development,” 

Part 2, Technical Proposal No. 71 M-l 173/C4163, Malibu, Calif., March 1971. 

2. Reed, George B., “Application of Kinematical Geodesy for Determining the Short Wave 
Length Components of the Gravity Field by Satellite Gradiometry,” The Ohio State 
University Research Foundation, Report No. 201, March 1973. 


20 


NASA- Langley, 1976 



NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 
WASHINGTON, D.C. 20546 


OFFICIAL BUSINESS 

PENALTY FOR PRIVATE USE *300 SPECIAL FOURTH-CLASS RATE 

BOOK 


POSTAGE AND FEES PAID 
NATIONAL AERONAUTICS AND 
SPACE ADMINISTRATION 
451 



632 001 Cl U E 760709 
DEPT OF THE &IH FCBCE 
AF WEAPONS LABORATORY 

&ttn : technical libb&by 
KTRTLAND AFB NM 87117 


S00903DS 

(SDL) 


POSTMASTER : 


If Undeliverable (Section 158 
Postal Manual) Do Not Return 


"The aeronautical and space activities of the United States shall he 
conducted so as to contribute ... to the expansion of human knowl- 
edge of phenomena in the atmosphere and space. The Administration 
shall provide for the widest practicable and appropriate dissemination 
of information concerning its activities and the results thereof 

— National Aeronautics and Space Act of 1958 


NASA SCIENTIFIC AND TECHNICAL PUBLICATIONS 


TECHNICAL REPORTS: Scientific and 
technical information considered important, 
complete, and a lasting contribution to existing 
knowledge. 

TECHNICAL NOTES: Information less broad 
in scope but nevertheless of importance as a 
contribution to existing knowledge. 

TECHNICAL MEMORANDUMS: 
Information receiving limited distribution 
because of preliminary data, security classifica- 
tion, or other reasons. Also includes conference 
proceedings with either limited or unlimited 
distribution. 

CONTRACTOR REPORTS: Scientific and 
technical information generated under a NASA 
contract or grant and considered an important 
contribution to existing knowledge. 


TECHNICAL TRANSLATIONS: Information 
published in a foreign language considered 
to merit NASA distribution in English. 

SPECIAL PUBLICATIONS: Information 
derived from or of value to NASA activities. 
Publications include final reports of major 
projects, monographs, data compilations, 
handbooks, sourcebooks, and special 
bibliographies. 

TECHNOLOGY UTILIZATION 
PUBLICATIONS: Information on technology 
used by NASA that may be of particular 
interest in commercial and other non-aerospace 
applications. Publications include Tech Briefs, 
Technology Utilization Reports and 
Technology Surveys. 


Details on the availability of these publications may be obtained from: 

SCIENTIFIC AND TECHNICAL INFORMATION OFFICE 
NATIONAL AERONAUTICS AND SPACE ADMINISTRATION 

Washington, D.C. 20546 




